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ABSTRACT 

We measure the radii and two-dimensional light prohles of a large sample of 
young, massive star clusters in M83 using archival HST/WFC3 imaging of seven 
adjacent helds. We use GALFIT to £t the two-dimensional light profiles of the 
clusters, from which we find effective (half-light) radii, core radii, and slopes of 
the power-law (EFF) profile (77). We find lognormal distributions of effective 
radius and core radius, with medians of ^2.5 pc and ?«1.3 pc, respectively. Our 
results provide strong evidence for a characteristic size of young, massive clusters. 
The average effective radius and core radius increase somewhat with cluster age. 
Little to no change in effective radius is observed with increasing galactocentric 
distance, except perhaps for clusters younger than 100 Myr. We find a shallow 
correlation between effective radius and mass for the full cluster sample, but 
a stronger correlation is present for clusters 200-300 Myr in age. Finally, the 
majority of the clusters are best fit by an EFF model with index 77 < 3.0. There 
is no strong evidence for change in 77 with cluster age, mass, or galactocentric 
distance. Onr results suggest that clusters emerge from early evolution with 
similar radii and are not strongly affected by the tidal held of M83. Mass loss 
dne to stellar evolntion and/or GMC interactions appear to dominate cluster 
expansion in the age range we stndy. 

Subject headings: galaxies: individnal: M83, galaxies: star clusters: general 


1. Introduction 

Star clusters are an important end result of the star formation process. Because they 
are gravitationally-bound objects, they can potentially survive for a long time and act as 
tracers of past star formation events. Dne to their high densities and surface brightnesses, 
massive star clusters can also be detected ont to large distances. In fact, with the Hubble 
Space Telescope (HST), individnal stars in clusters can be resolved beyond the Local Gronp 
(~5 Mpc away, Larsen et ah 2011), and the extended light profiles of clusters can be resolved 
ont to a few tens of Mpc (e.g., Whitmore et ah 1999; Schweizer 2004; Bastian et ah 2013). 
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Star clusters can also be approximated as simple stellar populations, so their properties (e.g., 
age and mass) can be determined fairly well from integrated photometry (e.g., Cabrera-Ziri 
et al. 2014). A signihcant fraction of stars are known to form in clustered environments 
(Lada & Lada 2003), and some fraction of these end up in bound clusters, depending on the 
intensity of star formation (Goddard et al. 2010; Adamo et al. 2011, 2015). Understanding 
cluster evolution and dissolution, which are connected to their sizes and concentrations, is 
therefore important for nnderstanding star formation in general. 

Populations of yonng massive star clusters (YMCs) in nearby galaxies have been the 
focus of many studies in recent years. Particularly, there is much interest in nnderstanding 
whether these objects are yonng analogs of globular clusters (GCs). If so, we could study 
this early mode of star formation at much higher resolution locally, and better nnderstand 
how GGs relate to their host galaxies. One intriguing connection between YMGs and GGs is 
the apparent uniformity in size across a wide range of age, mass, and environment. Stndies 
of star clusters in the Milky Way (e.g., Harris 1996; Portegies Zwart et al. 2010), Local 
Gronp (e.g., Elson et al. 1987; Mackey & Gilmore 2003b; Barmby et al. 2007), and other 
nearby galaxies (e.g., Larsen 2004; Jordan et al. 2005; Scheepmaker et al. 2007) have found 
the typical effective radius, or half-light radius, to be 2-3 pc. The lack of a signihcant 
relationship between cluster mass and radius (Zepf et al. 1999; Larsen 2004; Scheepmaker 
et al. 2007) is especially pnzzling, given that giant molecular clouds (GMGs) do have such 
a relation (Larson 1981). In addition, cluster age is apparently not strongly correlated with 
radius (e.g., Larsen 2004), and the hndings are mixed for galactocentric distance versus 
radius (e.g., van den Bergh et al. 1991; Jordan et al. 2005). The core radii of clusters can 
also be studied out to distances of a few Mpc (Schweizer 2004), and the spread of cluster 
core radii has been found to increase signihcantly with increasing age (Elson 1991; Mackey & 
Gilmore 2003b; Glatt et al. 2009). Finally, the power-law slope of young cluster light prohles 
have been measured in nearby galaxies, and some evidence of steepening or trnncation in 
the outer prohle wings with increasing age has been found (Whitmore et al. 1999; Larsen 
2004; Schweizer 2004). 

M83 is a face-on spiral galaxy [i = 27°, from the major to minor isophotal diameter 
ratio given in de Vaucouleurs et al. 1991) located at 4.5 Mpc (Thim et al. 2003). Several 
recent studies have scrutinized many aspects of M83’s substantial young cluster population 
(Ghandar et al. 2010; Bastian et al. 2011, 2012; Silva-Villa et al. 2013, 2014; Ghandar et al. 
2014; Hollyhead et al. 2015; Adamo et al. 2015). In particular, Bastian et al. (2012) measured 
the effective radii of a sample of young clusters located in two adjacent HST helds. The 
authors performed a careful visual inspection of their clusters to ensure that the majority 
of their objects are older than a crossing time, which implies that they are bound systems. 
The distribution of radii is lognormal, peaking at ~2.5 pc, and a slight expansion is seen 
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with age. 

In this paper, we build upon Bastian et ah (2012) by studying the sizes and light prohles 
YMCs in M83 in greater detail. We carefully select massive, isolated clusters with well- 
behaved light prohles from seven adjacent HST helds to ht with GALFIT (Peng et ah 2002, 
2010). This sample is one of the most extensive from a single galaxy for which the slope of 
the light prohle was left as a free parameter rather than set to a hxed value. We compare the 
effective radii, core radii, and slopes of the power-law light prohles of our clusters to their 
ages, masses, and galactocentric distances to probe the physical processes that dominate 
their structural evolution. We compare the results from our sample to other studies of YMC 
and GC sizes in the Milky Way, Local Group, and other nearby galaxies. 

This paper is organized as follows. We describe our HST observations and star cluster 
selection criteria in Section 2. In Section 3, we describe how we have measured the sizes 
and light prohles of the star clusters in our sample. In Section 4, we present the ehective 
radii, core radii, and light prohle slopes of the clusters and how they compare to other 
cluster properties. We discuss the implications of these results in Section 5, and hnally, we 
summarize this study in Section 6. 


2. Observations Cluster Catalogue 

We use archival HST/WFC3 imaging data retrieved from the Mikulski Archive for 
Space Telescopes (MAST), which were calibrated and drizzled on-the-hy when requested for 
download. These data were obtained by two programs, 11360 (PI; O’Connell) and 12513 
(PI; Blair). In total, seven helds were observed across the disk of M83 between August 2009 
and September 2012. See Figure 1 in Silva-Villa et ah (2014) for the location of the HST 
helds superimposed on M83. All seven helds were imaged in the F336W, F438W, F547M, 
and F814W hlters, except Field 1, for which F555W replaced F547M. 

The cluster catalogue from which our sample is selected is presented in Silva-Villa et ah 
(2014), further discussed in Adamo et ah (2015), and briehy summarized here. First, sources 
were selected from the F547M band in Fields 2 through 7 and F555W in Field 1 using Source 
Extractor (Bertin & Arnouts 1996). A concentration index (Cl) cut in F547M/F555W was 
applied to remove unresolved sources. The Cl is the magnitude diherence between aperture 
radii of 1 and 3 pixels, and it provides a measure of the relative size of an object. The 
Cl of an extended object is larger than that of an individual star, so applying a Cl cut 
is one method for removing stars from a star cluster catalogue (e.g., Holtzman et al. 1996; 
Whitmore et al. 2010). Sources located within 500 pc of the center of M83 were excluded due 
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to incompleteness because of crowding and regions of high extinction. Each resolved source 
was visually inspected and classihed as either a bona hde cluster (class 1), group/association 
(class 2), or false detection (class 3). Aperture photometry in the remaining hlters was 
performed on the catalogue. For each source with photometric detections and moderate 
errors (<0.3 mag) in all four hlters, the SED-htting technique presented in Adamo et ah 
(2010) was used to determine the age, mass, and extinction. 

In this study, we use the F547M/F555W images for htting the cluster light prohles. 
These hlters are not as strongly ahected by dust as compared to the bluer hlters, nor by 
bright (red) stars as compared to F814W. We limit our cluster sample to those sources in 
the Silva-Villa et al. (2014) catalogue of mass >10^ Mq and visual inspection class 1 (i.e., 
the most centrally concentrated and isolated cluster candidates; i.e., bona hde clusters). The 
mass limit minimizes the ehects of stochastic sampling of the stellar IMF, which can strongly 
ahect the ages and masses derived from broad-band magnitudes of low mass clusters, but 
diminishes at high masses (e.g., Popescu & Hanson 2010; Fouesneau & Langon 2010; Silva- 
Villa & Larsen 2011). Limiting our sample to clusters of class 1 ensures that their light 
prohles are well-resolved and reasonably well-behaved, which is optimal for measuring radii. 

Finally, a number of tests showed that clusters without extraneous sources within the 
htting region, such as bright stars or background huctuations, were the clusters that returned 
the most reliable, stable hts^. We therefore perform an additional visual inspection step to 
determine the degree to which a cluster is isolated within the 30 x 30 pixel htting region. 
To do this, we displayed circles of radius 15 pixels centered on each cluster that satished 
the mass and class limits discussed previously. An isolation hag, called /iso, was assigned to 
each cluster based on the number and brightness of sources within the circle. The hags are 
dehned as follows: 

1. Isolated. No nearby sources, or sources are faint enough to be indistinguishable from 
sky background. 

2 . Somewhat isolated. One or a few faint sources, or slight huctuations in sky background. 

3. One or two bright sources nearby. 

4. Crowded. Surrounded by several bright sources. 

In the analysis that follows, we report GALFIT results for only those clusters with /iso = 1 or 
2. Limiting our sample to relatively isolated clusters removes most of the youngest clusters 


^Also discussed on the GALFIT website: http://users.obs.carnegiescience.edU/peng/work/galfit/TFAQ.html#sensitivity 
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(less than a few tens of Myr). Such clusters tend to still be associated with their natal star- 
forming regions, and are therefore surrounded by young, bright stars that may be unrelated 
to the clusters themselves (Larsen 2004). The number of clusters remaining in the sample 
after the mass cut, visual inspection class cut, and isolation flag cut is 509. 

Our cluster sample is undoubtedly incomplete. It was selected using complicated meth¬ 
ods from images that cover a limited region of M83 and vary in exposure time. The clusters 
we study range in age, mass, surface brightness, and are subject to different local envi¬ 
ronments (e.g., crowding, dust lanes), though we have attempted to select only the best 
candidates for this study. For these reasons, the completeness of our sample is nearly impos¬ 
sible to quantify, but we indicate the ways in which we expect our sample to be incomplete 
in the following sections. 

A small amount of processing was required to prepare the HST images for use with 
GALFIT. First, GALFIT prefers input images to be in units of counts or electrons instead 
of flux units. We multiplied the drizzled images by the exposure time in order to convert 
from electrons per second to electrons. Because GALFIT multiplies by the GAIN keyword 
value in the image header to convert an input image from counts to electrons, and our images 
were already in units of electrons, we set the GAIN keyword equal to 1.0 e“/ADU. We also 
updated the image headers to include the approximate WFG3/UVIS readnoise, 3.11 e“, in 
the header keyword RDNOISE. 


3. Measuring Star Cluster Sizes 
3.1. Methodology 

We use the two-dimensional fitting package GALFIT (Peng et ah 2002, 2010) to fit 
the light prohles of the star clusters. Since the clusters are relatively young, we assume 
that an EFF prohle (also known as a Moffat prohle) accurately describes their true light 
prohles (Elson et ah 1987; Larsen 1999; Whitmore et ah 1999; Mackey & Gilmore 2003b). A 
comparison between hts to EFF, King (King 1966), and Wilson (Wilson 1975) prohles was 
performed by McLaughlin & van der Marel (2005) for YMGs and GGs in the Milky Way, 
Fornax dwarf galaxy, and Magellanic Glouds. They found that the more extended prohles, 
EFF and Wilson, ht young and old clusters better than the King model, in most cases. For 
most young clusters, the EFF prohle hts about as well as the Wilson prohle, though in some 
cases it does worse. We choose to go with the simpler EFF prohle because the Wilson prohle 
is not included in GALFIT, and McLaughlin & van der Marel (2005) also hnd that most 
of the physical properties of most clusters are well-constrained no matter the choice of light 
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profile. 

The EFF profile takes the form 


/i(r) = /io(l +rVa^) 


( 1 ) 


where /i is the surface brightness in a given band, a is a characteristic radius, and rj is the 
power-law exponent that determines the steepness of the prohle wings. Note that rj here is 
equal to 7/2 in Equation 1 of Elson et al. (1987). The core radius of a circular EFF prohle, 
which is the radius at which the surface brightness is half its central value, is given by 

reore = FWHM /2 = 2aV2i/r) _ 1. (2) 


The effective radius (refr), or half-light radius, is dehned to be the radius of the circular area 
which contains half of the total surface brightness of the light prohle. This can be written 
as _ 


reff = 2rc 


( 1 / 2 )^ 


- 1 


( 3 ) 


For an elliptical prohle, the true core and ehective radius can be found by multiplying 
Equations 2 and 3 by a factor of 0.5(1 -|- b/a), where h/a is the semiminor to semimajor axis 
ratio (ishape manual, Larsen 1999). 


For each cluster, we ht two components, an EFF component and a sky background 
component. Both are ht simultaneously over a 30 x 30 pixel box {^26 x 26 pc; 1 WFC3/UVIS 
pixel ~0.87 pc at the distance of M83). We assume the sky background is hat, but leave the 
amplitude of the sky as a free parameter. For the EFF component, all of the parameters 
are left free. These parameters are the x and y image coordinates of the cluster center, 
total magnitude, FWHM, 7 , axis ratio {b/a), and position angle. Table 1 details the initial 
guesses for each of the free parameters, excluding the x and y image coordinates which come 
directly from the Silva-Villa et al. (2014) catalogue. GALFIT returns the best-ht values for 
each free parameter and their la uncertainties. GALFIT returned an error for 31 out of the 
509 clusters in our sample, so these objects were excluded from the following analysis. 


In order to extract structural components from images, GALFIT convolves a model 
image with a PSF and compares the result to the data. An accurate PSF is therefore 
essential for reproducing the ehects of the telescope optics in the model images. We create 
PSFs for each held and hlter combination from several bright, isolated stars in each image 
by using pstselect and psf within DAOPHOT in IRAF. We spatially subsample these 
empirical PSFs by a factor of 10. 



3.2. Testing the Initial Gnesses 


We investigate the effect of varying the input parameters by performing tests on the most 
isolated, well-behaved clusters in Field 5 (/iso = !)• We varied one of the following initial 
guesses at a time: FWHM, 77 , and size of the htting box. The remaining input parameters 
were kept the same as in Table 1. For each cluster, we calculated the difference between 
the best-£t FWHM and rj from the GALFIT run that used the standard initial guesses 
and the runs in which we varied a single parameter. We denote these differences AFWHM 
and Af]. The average and standard deviation of these differences (e.g., (AFWHM) and 
(tAFWHM) are listed in Table 2 for each varied input parameter. We find that changing the 
initial guesses and the htting box size have little affect on the best-£t FWHM. The average 
FWHM differences are close to zero, and the scatter about the average is small. The average 
differences of the best-£t rj are also close to zero for all initial guesses an box sizes, but the 
scatter tends to be larger than that of the best-£t FWHM. Changing the htting box size also 
results in average diherences close to zero, but there is a larger scatter about those averages, 
especially for rj. Therefore, although changing the htting box size or initial guesses may 
ahect hts to individual clusters, on the whole, the best-ht results are essentially the same. 


4. Results 

In this section, we compare the structural parameters from GALFIT with other cluster 
properties. We also present these values in Table 3 for the 478 clusters that were well-ht by 
GALFIT. As discussed in Sections 4.1 and 4.2, the ehective radii of clusters best-described 
by 77 < 1.3 are not well-constrained, and all structural parameters for clusters best-ht by 
77 < 1.1 should be treated with caution. 


Table 1: GALFIT Input Parameters 


Parameter 

Value 

Total Magnitude 

20.0 mag 

FWHM 

2.5 pix 

V 

1.5 

Axis Ratio 

1.0 

Position Angle 

25° 

Sky Background 

300.0 e- 
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Table 2: Tests of GALFIT Input Parameters 


Input 

FWHM (pix) 

^clusters 

(AFWHM) (pix) 

Output 

aAFWHM (pix) 

(At/) 

aAr] 

1.5 

86 

-0.007 

0.08 

0.01 

0.4 

3.5 

86 

-0.01 

0.09 

0.006 

0.4 

V 



1.1 

86 

0.03 

0.1 

0.1 

0.5 

2.0 

86 

-0.005 

0.08 

-0.04 

0.2 

Fitting Box (pix) 



20 x20 

84 

0.04 

0.3 

0.05 

0.8 

40x40 

84 

-0.01 

0.2 

0.003 

0.6 


4.1. Effective Radius 

Measuring is difficult when the power-law index of the EFF profile, rj, is near 1.0. 
When ?7 < 1.0, rgfj is undefined (unless an outer truncation radius is assumed, as in Larsen 
et ah (2001) and Larsen (2004)). When rj is slightly greater than 1.0, rgg diverges to un¬ 
physical values, as does the error, Sves, as shown in Figure 1. In this plot, we show the 
fractional error of refr with respect to rj, assuming 10% errors in FWHM, rj, and axis ratio 
h/a. In order to avoid unphysical effective radii and errors, we study clusters that are best 
fit by 7] > 1.3, which corresponds to a fractional error in effective radius of ~50%, in this 
section. We briefly discuss the biases this r] limit imposes on our sample in Section 4.1.1. 
The number of clusters remaining after including this r] limit is 247, which is about half of 
the 478 clusters that were well-fit by GALFIT. 

In addition, we also impose an age limit on our sample of <300 Myr because the cluster 
age estimates from SED-fltting become fairly uncertain above this age. Stellar evolution¬ 
ary models that use different AGB star prescriptions show significant discrepancies around 
300 Myr. Therefore, the cluster age determination of the SED fits becomes strongly model 
dependent at this point. In addition, the age-metallicity degeneracy of globular clusters 
(GGs) comes into play. Metal-poor GGs a few Gyr old can populate an area in color space 
that is also occupied by 300 Myr solar-metallicity clusters. The number of clusters with ages 
<300 Myr is 152, and the number with ages >300 Myr is 95. 



Table 3. Properties of YMCs in M83 
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Fig. 1.— Fractional error on the effective radius as a function of the power-law index, rj. To 
calculate Sr^s, we assume 10% errors on FWHM, t], and the axis ratio h/a. As shown by the 
dashed line, the fractional error is ~50% at ?7 ~ 1.3. 

% i. i. Distribution of Effective Radii 

In Figure 2, we present the distribution of effective radii for the conservative rj limit 
discussed above, rj > 1.3 (top panel), and a less conservative limit, ?7 > 1.1 (bottom panel). 
Clusters younger than 300 Myr are plotted in blue and clusters older than 300 Myr are plotted 
in gray. In the top panel, both age groups are well-described by lognormal distributions, as 
shown by fits to the data plotted as solid and dashed black lines. The median effective radius 
of the young sample in the top panel is ^2.5 pc, while that for the old sample is Ri2.8 pc. 
The full range of effective radii covered by either sample is ~0.3 to 10 pc. 

Including clusters best-fit by smaller r] values introduces a tail of more extended objects 
into the distribution, as seen in the bottom panel. The presence (and location) of the peak 
in the distribution does not change when these extended objects are introduced. In addition, 
only a few very compact clusters (<0.4 pc) appear in either panel due to the Cl cut described 
in Section 2. We can easily distinguish clusters 1 to 1.5 pc in size from unresolved objects, so 
the dropoff from the peak to smaller sizes in either panel is not due to biases in the sample. 
The median effective radius is ~3.0 pc for the young sample and ~3.5 pc for the old sample 
in the bottom panel. We therefore find strong evidence for a characteristic size of YMCs in 
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log(reff/pc) 

Fig. 2.— Distribution of effective radii. Clusters younger than 300 Myr are plotted in blue 
and clusters older than 300 Myr are plotted in gray. Top panel: Only clusters best-£t by 
> 1.3 are included. Lognormal fits to the two age groups are plotted as solid (young 
clusters) and dashed (old clusters) black lines. Bottom panel; Clusters best-fit by 77 > 1.1. 
The error bars are median errors of clusters located in 0.5 dex bins in effective radius. 

M83 of about 2.5 to 3 pc. 

We also plot median error bars for each 0.5 dex in effective radius in the bottom panel to 
show that the extended clusters have much larger errors than the smaller clusters. Therefore, 
while there does exist a population of extended clusters in M83, the errors are too large to 
include them in our analysis. In the remainder of Section 4.1, we only consider clusters 
described by rj > 1.3. 

The distribution in the top panel of Figure 2 is very similar to that found by Bastian 
et ah (2012) in Fields 1 and 2 of M83 (their Figure 13). This includes the overall lognormal 
shape, peak radius, and range of radii present. These similarities are clear despite the fact 
that Bastian et ah (2012) used ISHAPE (Larsen 1999) instead of GALFIT, and assumed 
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a fixed rj of 1.5. It seems that allowing rj to remain a free parameter does not drastically 
effect the shape of the distribution. In other words, the effective radius appears to be fairly 
accurately recovered from a well-behaved light prohle by a two-dimensional £t, even if the 
assumed shape of the light prohle is somewhat inaccurate, as previously suggested by Kundu 
& Whitmore (1998), Larsen (1999), and McLaughlin & van der Marel (2005). 

The distribution in Figure 2 is also quite similar to that found by studies of young 
massive clusters (Meurer et ah 1995; Larsen 1999; Whitmore et ah 1999; Larsen 2004; Bastian 
et al. 2005; Lee et al. 2005; Barmby et al. 2006; Scheepmaker et al. 2007; Mayya et ah 2008; 
Portegies Zwart et al. 2010; Annibali et al. 2011; Johnson et al. 2012; San Roman et al. 2012) 
and globular clusters (Jordan et al. 2005; Spitler et al. 2006; Georgiev et al. 2008; Harris 
2009; Harris et al. 2010; Masters et al. 2010) in the Milky Way and other nearby galaxies. 
Some galaxies also harbor populations of extended clusters, such as the “faint fuzzies” in 
NGC 1023 (Larsen & Brodie 2000), the extended globular clusters in M31 (Huxor et al. 2005), 
and the extended, faint clusters in MlOl (Simanton et al. 2015). The effective radii of these 
objects range from ~7 to ~30 pc, but they may not be all that different structurally than 
extended GGs in the Milky Way’s outer halo (van den Bergh & Mackey 2004). The “faint 
fuzzies”, however, appear to be a unique population of objects, which have very different 
kinematics, metallicities, and luminosities than extended Milky Way GGs (Brodie & Larsen 
2002). Massive clusters of all ages in the Magellanic Glouds also appear to be somewhat 
extended compared to those in most other systems, having median effective radii of ~8 pc 
(McLaughlin & van der Marel 2005; Glatt et al. 2009). We discuss these populations of 
extended clusters further in Section 5.2. 


4-1-2. Effective Radius versus Age 

In Figure 3, we plot effective radius as a function of age. For bins spanning 0.3 dex in 
age containing more than 5 clusters, we calculate the median and mean effective radii, which 
are plotted as orange triangles and yellow squares, respectively. We also plot one standard 
deviation about the mean for each bin as dashed orange curves. Finally, we perform a least- 
squares £t to the mean effective radii, including the standard errors on the mean, and plot 
the result as a black line. We find a clear correlation between mean radius and age with 
a slope of 0.26 ± 0.07. The mean mass of the clusters in each bin is relatively constant at 
~20,000 Mq, so this correlation is unrelated to cluster mass. 

As we discuss in Section 2, our cluster sample suffers from incompleteness which is 
difficult to quantify. The most obvious bias in our sample is the lack of young clusters, <30- 
50 Myr in age, due to the isolation hag cut discussed in Section 2. Of the young clusters 
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Fig. 3.— Effective radius as a function of cluster age. The blue points are clusters <300 Myr 
in age and best-fit by r; > 1.3. Median and mean effective radii are calculated in 0.3 dex age 
bins containing 5 or more clusters, and are plotted as orange triangles and yellow squares, 
respectively. One standard deviation in effective radius about the mean is plotted as dashed 
orange curves. The black line is a least-squares fit to the mean effective radii, which includes 
the standard errors on the mean. 

that exist in our sample, we would expect that all sizes are well-represented because young 
clusters are relatively easy to detect due to their brightness. On the other hand, it is likely 
that extended, old clusters are not fully represented. The low surface brightness of these 
objects make them difficult to detect against the background. Thus, the trend in Figure 3 
could actually be steeper if there is a signihcant population of faint, extended clusters that 
has gone undetected. 

This trend is similar to the size-age correlation found by Bastian et al. (2012), though 
our data show larger scatter and cover a smaller range in age. In fact, Bastian et al. (2012) 
note that the correlation they find may be due to the fixed power-law index of the light 
prohle, but that hypothesis is now ruled out with this study. Several other studies have 
also looked for a relation between effective radius and age in young massive clusters. Larsen 
(2004) Ends essentially no trend for a large sample of clusters from several nearby galaxies. 
Barmby et al. (2009) also End no trend for 23 clusters in M31, and San Roman et al. (2012) 
End no trend for 161 clusters in M33. The young clusters studied by McLaughlin & van der 
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Marel (2005) in the LMC and SMC show no relationship between age and effective radius. 
In M51, three studies hnd a slight positive relationship with a power-law slope of ~0.1 (Lee 
et ah 2005; Scheepmaker et al. 2007; Hwang & Lee 2010). On the other hand, another study 
finds a very slight negative trend in M51, though clusters smaller than 2 pc were unresolved 
and not included in the £t (Bastian et al. 2005). 


4.1.3. Effective Radius versus Galactocentric Distance 

In Figure 4, we plot effective radius as a function of galactocentric distance. Again, 
we only include clusters <300 Myr in age and best-fit hy t] > 1.3. The least-squares fit 
to the mean effective radius in each 0.2 dex bin hnds a correlation with slope 0.29 ± 0.09. 
However, cluster age is also weakly correlated with galactocentric distance, such that within 
2 kpc of the center of M83, clusters are ^140 Myr old on average. The average age increases 
to ?^230 Myr beyond 4 kpc from the center. We therefore separate the clusters into three 
100 Myr bins in age in Figure 5 to remove the age-distance relation. We perform least- 
squares fits to the mean effective radii in each age bin, and plot the results as black lines 
with the slopes labeled in the lower right hand corner of each panel. The bottom two panels 
show little to no correlation between mean radius and galactocentric distance, while the top 
panel, containing the youngest clusters, shows a moderate correlation. However, there is 
considerable scatter in this plot, so it is unclear if the trend is real. 

Past studies of young massive clusters have found mixed results with this relationship. 
Bastian et al. (2012) notes that clusters in Field 2 (further from the center) are on average 
0.4 pc larger than those in Field 1 (closer to the center). They attribute this difference to the 
relationship between cluster age and radius, which corresponds with our results. In M51, the 
evidence is mixed: Bastian et al. (2005) find no obvious trend between projected galactocen¬ 
tric distance and effective radius, but Scheepmaker et al. (2007) find a weak relation with a 
power-law slope of 0.12 ±0.02. In M33, there is no correlation between galactocentric radius 
normalized by circular velocity and cluster size, unless a small number of large, outer clusters 
are included (San Roman et al. 2012). There is a somewhat stronger trend in NGC 7252, 
where Bastian et al. (2013) find a slope of 0.35 ± 0.20 for very massive clusters, >10® M©. 

For GCs, the picture is almost the same, though most studies cover a much larger range 
in galactocentric distance than our sample. In the Milky Way, there is a strong relationship 
between galactocentric distance and radius, van den Bergh et al. (1991) find a power-law 
slope of 0.5. In contrast, M31 shows little to no correlation (Barmby et al. 2007). Studies 
of GGs in nearby early-type galaxies find much shallower relations (~0.1), or almost none 
at all (Jordan et al. 2005; Madrid et al. 2009; Harris 2009; Harris et al. 2010; Masters et al. 
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Fig. 4.— Effective radius as a function of galactocentric distance. Points and lines are the 
same as in Figure 3, except that the mean, median, and standard deviation are found for 
0.2 dex bins in galactocentric distance. 

2010), although Spitler et ah (2006) hnd a somewhat stronger relation for a GCs around the 
Sombrero Galaxy. Further, in dwarf irregular galaxies, Georgiev et al. (2008) hnd no clear 
trend between galactocentric distance and cluster size. 


4-1.4- Effective Radius versus Mass 

As discussed in Gieles et al. (2006), disruption of YMGs by molecular clouds depends on 
cluster density. Quantifying the mass-radius relation is therefore important for understand¬ 
ing cluster disruption in gas-rich environments like galaxy disks (e.g.. Fall & Ghandar 2012). 
In Figure 6, we plot effective radius as a function of cluster mass. We perform a least-squares 
£t to the mean effective radii in bins of 0.2 dex in mass and hnd a shallow correlation with 
slope 0.3 ± 0.1. To remove the known age-radius relation shown in Section 4.1.2, we divide 
the young clusters into three 100 Myr bins in age and perform a least-squares ht to the 
clusters in each bin, as shown in Figure 7. The correlation between mass and mean ehective 
radius becomes stronger with increasing cluster age. Notably, the 200-300 Myr bin has a 
slope of 0.5 ± 0.2, although the range of cluster masses contained in this bin is fairly narrow 
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Fig. 5.— Same as Figure 4, except that the cluster sample is now divided into three bins in 
age. 

and may be skewing the results. 

As discussed previously, faint, extended clusters are likely missing from our sample 
because they are particularly difficult to detect. However, the correlation between mean 
radius and mass we hnd for clusters 100-300 Myr in age would not necessarily be diminished 
unless a large population of low-mass, extended clusters has gone undetected. We expect the 
sample to be dehcient in extended clusters of all masses, not necessarily low-mass clusters 
alone. 

Many studies of young massive clusters have found little to no correlation between 
effective radius and mass (or luminosity) (Meurer et ah 1995; Zepf et ah 1999; Larsen 2004; 
Lee et al. 2005; Bastian et ah 2005; Scheepmaker et ah 2007; Mayya et ah 2008; Barmby 
et al. 2009; Hwang & Lee 2010; Bastian et al. 2012). There are a few interesting exceptions. 
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Fig. 6.— Effective radius as a function of cluster mass. Points and lines are the same as in 
Figure 3, except that the mean, median, and standard deviation are found for 0.2 dex bins 
in mass. 

however. The data provided by McLaughlin & van der Marel (2005) for intermediate-age 
clusters (>300 Myr) and GCs in the Magellanic Clouds show a rather strong relation, on the 
order oc (see Figure 4 in Fall & Chandar 2012). Bastian et ah (2013) find a power- 
law slope relation of 0.29 ±0.09 for the very massive (>10® Mq) young clusters in NGC 7252, 
and little to no relation below this mass. Kissler-Patig et ah (2006) hnd a similar relation for 
a larger sample of very massive young clusters. GCs show little to no relationship between 
luminosity and radius (van den Bergh et ah 1991; Kundu & Whitmore 2001; Madrid et ah 
2009), except for perhaps the most luminous (most massive) objects, which may have larger 
radii (Spitler et al. 2006; Barmby et ah 2007; Harris 2009; Harris et al. 2010). Therefore, for 
all types of clusters, there is essentially no evidence for a strong mass-radius relation below 
~10® Mq, but higher mass objects do appear to get larger with increasing mass. 


4.2. Core Radius 

Considering the core radius is advantageous because it is less biased against extended 
clusters than the effective radius. The best-fit FWHM from GALFIT is equal to twice the 
core radius, as described in Section 3.1. Because the core radius does not depend on r^, unlike 
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Fig. 7.— Same as Figure 6 , except that the cluster sample is now divided into three bins in 
age. 

the effective radius, we do not have to impose the same r] limit as we do in Section 4.1. We 
relax the limit to 77 > 1.1 in this section for the following reason. GALFIT cannot fit the light 
profiles of clusters with r] < 1.0, which by definition contain an infinite amount of light. Even 
though this is unphysical, Larsen (2004) has shown that clusters with < 1.0 light profiles 
do exist. We find that GALFIT assigns these ‘unphysical’ clusters r] values slightly higher 
than 1.0. This r] limit does remove some clusters that are accurately described by 77 ~ 1.0, 
but since GALFIT cannot distinguish these clusters from those with unphysically shallow 
profiles, we remove all clusters best-fit by 77 < 1.1. With this somewhat relaxed restriction, 
the number of clusters <300 Myr in age increases to 224, and the number >300 Myr in age 
is 133. 

We plot the distribution of core radii in Figure 8 . Both age groups are fairly well- 
described by a lognormal distribution (with some skew to the left), as shown by fits to the 
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Fig. 8.— Distribution of core radii. Clusters younger than 300 Myr are plotted in blue 
and clusters older than 300 Myr are plotted in gray. Only clusters best-£t by r; > 1.1 are 
included. Lognormal fits to the two age groups are plotted as solid (young clusters) and 
dashed (old clusters) black lines. 

data represented by solid and dashed black lines. The median core radius of the young 
sample is f»1.3 pc and that for the old sample is f»1.6 pc. The full range of core radii is 
~0.1 to 4 pc. Even though this distribution is less biased against extended clusters, it is still 
strongly peaked, which is further evidence for a characteristic size of YMCs in M83. 

The distribution in Figure 8 is somewhat different than that found by Larsen (2004). 
In that study, the distribution of FWHM (equivalent twice our core radius) peaks at smaller 
values and declines more gradually to larger values than our distribution. We do, however. 
End a similar overall range of core radii. The core radii of 23 clusters in M31 span about the 
same range of values as Figure 8 as well (Barmby et ah 2009). Clusters in the Magellanic 
Clouds, however, seem to have larger average core radii than those found in nearby spirals 
(Elson et ah 1987; Elson 1991; Mackey & Gilmore 2003a,b; Glatt et ah 2009). 

A significant increase in the spread of core radii with increasing cluster age was observed 
for clusters in the LMG by Elson (1991) and Mackey & Gilmore (2003b), in particular. San 
Roman et ah (2012) End a similar result, though they measure smaller core radii on average, 
in M33. In Figure 9, we plot the core radii of our sample as a function of age (blue points). 
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Fig. 9.— Core radius as a function of age. The blue points are M83 clusters <300 Myr in 
age and best-fit hj rj > 1.1. The gray diamonds are LMC clusters from Mackey & Gilmore 
(2003b) of similar age to the M83 clusters. Median and mean core radii of the M83 clusters 
are calculated in 0.3 dex age bins containing 5 or more clusters, and are plotted as orange 
triangles and yellow squares, respectively. One standard deviation in core radius about the 
mean is plotted as dashed orange curves. The black line is a least-squares £t to the mean 
core radii of the M83 clusters, which includes the standard errors on the mean. 


and overplot LMC clusters of similar age from Mackey & Gilmore (2003b) for comparison 
(gray diamonds). Clearly, the core radii of the LMC clusters in this age range are larger on 
average than the clusters in M83. We do not see a strong increase in the spread of core radii 
in M83, but we do see that the core radius of the largest cluster increases with age, similar 
to the LMC. We hnd evidence for an overall trend of increasing core radius with increasing 
age. This trend has a slope of 0.24 ± 0.05, similar to the size-age trend in Figure 3. We 
do not see a bifurcation in Figure 9, which Mackey & Gilmore (2003b) suggest is present in 
their data. 
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Fig. 10.— Distribution of EFF power-law indices, rj. Clusters younger than 300 Myr are 
plotted in blue and clusters older than 300 Myr are plotted in gray. Only clusters best-£t 
by h ^ 1-1 cire included. 


4.3. EFF Index 

We now turn to the EFF prohle power-law index, p. Again, we have imposed an rj limit 
of ?7 > 1.1 for the reasons described in Section 4.2. In Figure 10, we plot the distribution 
of best-£t 7] values for our cluster sample. The distribution peaks at the lowest rj values, 
followed by a rapid falloff to higher r] values. The median of the young distribution is 1.5, 
similar to that found by previous studies (Elson et ah 1987; Mackey Sz Gilmore 2003b; Larsen 
2004; McLaughlin & van der Marel 2005; Glatt et ah 2009). A few objects are best-fit by 
rj > 3.0, up to 7] ~ 8.7, but we choose to focus on the region between 1.0 and 3.0 where most 
of our clusters lie. This range is consistent with the previously reported best-£t rj values 
for young clusters (Elson et ah 1987; Mackey & Gilmore 2003b; Larsen 2004; McLaughlin & 
van der Marel 2005; Glatt et ah 2009). Therefore, our choice to only £t relatively isolated 
clusters does not seem to bias our determination of the distribution of EFF power-law slopes. 

We also look for any change in rj with other cluster properties. In Figure 11, we divide 
the 7] distribution into three bins in cluster age in the left column, four bins in galactocentric 
distance in the middle column, and four bins in mass in the right column. We perform 
Anderson-Darling tests to determine if the distributions within each column are signihcantly 
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Fig. 11.— Distribution of EFF power-law indices divided into three bins in age (left column), 
four bins in galactocentric distance (middle column), and four bins in mass (right column). 
Only clusters <300 Myr in age and best-£t hj t] > 1.1 are included. 

different from one another. This takes the form of rejecting or accepting the null hypothesis, 
i.e., the distributions in each column could be drawn from the same parent distribution. For 
instance, for the cluster age column, we perform Anderson-Darling tests comparing the hrst 
bin to the second, the second to the third, the hrst to the third, and all three together. 
The same is done for the other two columns, though the number of permutations increases 
when a fourth bin is included in the tests. We report the significance level with which the 
null hypothesis can be rejected for a given test. The p-values of each test that rejects the 
null hypothesis are less than or equal to the significance level. The Anderson-Darling test is 
ideal for comparing the distributions in Figure 11 because it is a non-parametric test, which 
means it is not necessary to know the shape of the true underlying distribution. 
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For the bins in the age column, we cannot reject the null hypothesis for any bin as 
tested against any other bin. Therefore, the rj distributions in each bin are consistent with 
being drawn from the same parent population. We hnd no evolution of rj with age. 

For the bins in the galactocentric distance column, we can reject the null hypothesis at 
a significance level of 2% for each permutation that tests the distribution in the fourth bin 
(>4.0 kpc) against those in the three other bins. Therefore, the probability of obtaining as 
extreme of a difference as we hnd between the fourth bin and the other three bins is 2%. 
The top three bins are consistent with being drawn from the same parent population. We 
speculate that the reason for the distinctness of the distribution in the fourth bin is that it is 
sharply peaked at low values of rj, and contains very few clusters best-described by ?7 > 1.5, 
unlike the other three bins. It is also possible that this difference is driven by the lack of 
clusters detected beyond 4 kpc. 

For the bins in the mass column, we can reject the null hypothesis at a signihcance level 
of 4% for tests that compare the first bin (<15,000 Mq) to the third and fourth bins. Every 
other comparison between bins results in no rejection of the null hypothesis, i.e., the other 
combinations are consistent with being drawn from the same parent distribution. There 
are many fewer clusters above 20,000 Mq in our sample than clusters <15,000 M©, so this 
difference in sample size could be driving the test results. Therefore, while the Anderson- 
Darling tests suggest there may be signihcant differences between some of the distributions 
in Figure 11, the differences are rather small, and therefore should be further tested with 
studies of cluster populations in other galaxies. 

Larsen (2004) find a tentative increase of the mean rj (which they call a) with cluster 
age, although this is mostly due to their youngest age bin containing a larger fraction of 
clusters with 77 < 1.0 than the other age bins. Since we cannot determine which of our 
clusters are truly best described by 77 < 1.0, we cannot test this result. Using the data 
presented in Mackey & Gilmore (2003b) and McLaughlin & van der Marel (2005) for LMC 
and SMC clusters, we looked for a relationship between age and 77 and found none, although 
no clusters in those samples are consistent with 77 < 1.0. 


5. Discussion 

5.1. Drivers of Cluster Structural Evolution 

The expulsion of leftover gas from formation appears to be completed within the first 
4 Myr of a star cluster’s lifetime (Hollyhead et ah 2015). Excluding the two youngest 
(?^5 Myr) objects, the clusters in our sample are therefore likely to be gas-free. The radii we 
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measure for these clusters cannot be considered their initial radii, since clusters are expected 
to expand signihcantly between formation and the end of gas expulsion (e.g., Portegies Zwart 
et ah 2010). Similarly, the initial light profile slopes may have changed during this early phase 
of dynamical evolution. While we cannot know the initial conhgurations of our clusters, we 
can characterize their structural evolution between ~10 to 300 Myr. 

The dominant internal mechanism affecting the structure of clusters 10 to 300 Myr in age 
is mass loss due to stellar evolution (e.g., Portegies Zwart et ah 2010). If mass is lost slowly 
as compared to the dynamical time, an isolated cluster can maintain virial equilibrium while 
it gradually expands. Generally, mass loss can be assumed to be slow for massive clusters 
because even if a single star loses the majority of its mass very quickly, the amount of mass 
lost by the cluster is a relatively small fraction of the total (Reggie & Hut 2003). If the 
cluster is initially mass-segregated, this expansion will occur mostly in the core, since the 
massive stars in the core will lose relatively more mass than the rest of the population. On 
the other hand, the effective radius of a tidally-limited cluster should contract as it loses 
mass, but its core will expand (Heggie & Hut 2003). Therefore, we would expect to see an 
expansion of core radii, and an expansion or contraction of effective radius with age if mass 
loss through stellar evolution dominates over external forces. Since we hnd an expansion of 
the mean effective and core radii with cluster age in Figures 3 and 9, it seems likely that our 
clusters are not tidally-limited. 

Two-body relaxation also changes cluster structure, but typically becomes dominant 
over longer timescales than we encounter here (Portegies Zwart et al. 2010). We test this 
by calculating approximate half-mass relaxation times of the clusters using Equation 14.13 
from Heggie & Hut (2003). For a 100 Myr old, 10"^ Mq star cluster with an effective radius 
of 3 pc, the present-day, half-mass relaxation time is ^200 Myr. Here we assume that the 
three-dimensional half-mass radius is equal to 4/3 x the effective radius (Spitzer 1987). We 
also assume that the average stellar mass is 0.65 M©, which was estimated using a Kroupa 
IMF over a mass range of 0.1 to an approximate turn-off mass for a 100 Myr cluster, 6 M©. 
The relaxation time is strongly dependent on the cluster radius such that more compact 
clusters have shorter relaxation times. Given the distributions of radii and masses in our 
sample, the relaxation times range from a few tens of Myr to a few Gyr, with a peak around 
300 Myr. About 80% of our clusters are younger than one relaxation time, and the rest 
are about equal to or a few times older than their relaxation times. However, given that 
clusters are expected to expand drastically shortly after their formation, and that we observe 
a modest expansion during the hrst few hundred Myr of cluster evolution, it is likely that 
the half-mass relaxation times were shorter in the past. Whether two-body relaxation was a 
dominant process immediately after their formation, we cannot say, but it appears that the 
majority of the clusters in our sample are not strongly affected by two-body relaxation at 
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the present day. 

Two external processes that may influence cluster structure are the tidal held of M83 
and interactions with giant molecular clouds (GMCs). As mentioned previously, clusters that 
are tidally-limited behave differently during the period of significant stellar mass loss than 
those that are not. A cluster that is not tidally-limited can evolve essentially as an isolated 
system. The lack of a strong relationship between mean effective radius and galactocentric 
distance in Figures 4 and 5 suggests our sample is not strongly affected by tides. If it was, 
we would expect the clusters to become much larger with increasing galactocentric distance 
as the tidal held becomes weaker (e.g., Madrid et ah 2012). One potential signature of tidal 
influence on our sample is the tentative preference for shallower light profiles further from 
the galaxy center in Section 4.3. 


A further test of the influence of the tidal held on the clusters is to calculate approximate 
Jacobi radius, which defines the zero-velocity surface of a cluster in a tidal held. Stars within 
the Jacobi radius can be said to “belong” to the cluster. The Jacobi radius is defined by 
Equation 9 in Portegies Zwart et ah (2010), 


rj = 



( 4 ) 


where M is the mass of the cluster and u is the angular velocity of the cluster around the 
galactic center. We estimate u by assuming the angular velocity of molecular gas in the 
disk of M83 matches, to first-order, that of star clusters located at the same galactocentric 
distance. We found an approximate linear relation between uj of CO in M83’s disk and 
galactocentric distance from the top panel of Figure 4 in Lundgren et ah (2004). Using this 
relation, we estimate the Jacobi radius of each cluster in our sample, and And that they 
range in size from ~14 to ~43 pc. Above ~ 0.15 (Henon 1961) or ^0.2 (Alexander 

et ah 2014), a cluster Alls its Roche volume. We And Jacobi to half-mass radii ratios between 
~0.02 to ^0.5. About 68% of our clusters have r ^m /rT < 0.2, and 43% have rhm/''"j < 0.15. 
Therefore, based on rough estimates of the Jacobi radii, our sample is straddling the line 
between being influenced by tides and not being influenced by tides, i.e., some of our clusters 
are likely underfilling their Roche volumes while others are Ailing them. Comparing these 
results to cluster age, we And that larger Jacobi to half-mass radii ratios are more common 
in older clusters, so perhaps clusters become tidally-limited as they age. 


Cieles et al. (2006) studied the effect of interactions between CMCs and star clusters on 
the mass loss rate and disruption time of clusters. They show that a 10^ Mq cluster in the 
solar neighborhood would lose enough mass to be disrupted in about 2 Cyr. Each interaction 
leads to an energy gain by the cluster, though the magnitude of that gain is dependent on 
the characteristics of the cloud, the cluster, and the interaction itself. We would expect an 
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energy gain by the cluster to cause expansion (due to a shallower potential well), but it 
is unclear if the degree of tidal hlling would have an impact on this expansion. Assuming 
expansion occurs, and that it takes place at a similar rate to mass loss, the clusters in our 
sample may be signihcantly affected by GMC interactions over the age range we study. 


5.2. Similarity of Cluster Sizes 

In Section 4.1.1, we compare the radius distribution of our cluster sample with that of 
other samples in the literature. Our sample is one of the most extensive in a single galaxy 
for which the slope of the light prohle ( 77 ) was left as a free parameter, and yet we hud 
essentially the same distribution as found in many other studies. As many of these studies 
note, star clusters of a wide range in age, mass, and environment all seem to be ~3 pc in size. 
Even GCs, which have undergone a Hubble time of dynamical evolution and been affected 
by different physical processes (e.g., two-body relaxation, tidal shocks, core collapse), are 
similar in size to YMGs. 

In Section 4.1.1, we note the existence of populations of extended clusters (both young 
and old) in the Magellanic Glouds, M31, MlOl, several SO galaxies, and even the Milky 
Way, calling into question the conclusion that cluster sizes are relatively constant. In the 
Magellanic Glouds, clusters that have been carefully measured are, on average, larger than the 
clusters in M83 (see Figure 9). However, this well-dehned sample contains only 63 objects 
from a combined population of a few thousand in both Glouds (Glatt et ah 2010). The 
larger average size of Magellanic Gloud clusters could therefore be a selection effect, i.e., more 
compact clusters simply have not been studied at the same level of detail. Another possibility 
is that local environment somehow does play a role in determining cluster structure. Studies 
of larger samples of Magellanic Gloud clusters and of clusters in other dwarf galaxies are 
needed to answer this question. 

A number of extended clusters do exist in M83, as shown in the bottom panel of Figure 2, 
although we cannot precisely determine their sizes. We note, however, that the number of 
extended clusters is relatively small as compared to the rest of the population. This seems 
to be the case in most studies that hud extended clusters in other galaxies; clusters ~3 pc in 
size outnumber those of much larger size (e.g., Larsen & Brodie 2000; Simanton et al. 2015). 
These extended clusters appear to represent a tail in the size distribution rather than a 
separate population, and the majority of star clusters seem to have a characteristic compact 
size of ~3 pc. Faint fuzzies may be an exception, or could simply be old open clusters that 
are easier to detect in SO galaxies than star-forming spirals (Ghies-Santos et al. 2013). Of 
course, deeper surveys are needed to determine if large numbers of faint, extended clusters 



in nearby spiral galaxies have gone undetected. 

Perhaps the most surprising feature of the similarity between cluster sizes is the weak 
mass-radius relationship, although we find that this relationship strengthens as clusters age 
(Section 4.1.4). GMCs have a fairly strong mass-radius relationship, on the order of rcivic oc 
Mq^c (Larson 1981). It has been argued extensively that some process, either during cluster 
formation or subsequent evolution, must flatten the GMC mass-radius relationship. Ashman 
& Zepf (2001) suggest that if the star formation efficiency of GMGs depends on their binding 
energy, then the clusters that form from the GMGs will have no mass-radius relationship. 
This, however, implies that the subsequent dynamical evolution of the clusters cannot have a 
strong effect on their sizes. Gieles et ah (2010) explore the possibility that all self-gravitating 
stellar systems have the same mass-radius relation upon formation and lower mass clusters 
have moved away from this relation due to expansion driven by stellar mass loss and hard 
binaries. This idea is supported by studies briefly discussed in Section 4.1.4, which find a 
clear mass-radius correlation for objects >10® M©. In addition, the theory presented in Gieles 
et al. (2010) implies a similar rate of expansion as we observe in M83. However, the mass- 
radius relation we find becomes stronger with increasing cluster age, which is inconsistent 
with the suggestion by Gieles et al. (2010) that the initial mass-radius relation should degrade 
over time as low-mass clusters expand. 


6. Conclusions 

We have measured the effective radii, core radii, and EFF power-law index of a large 
sample of YMGs in M83. To do this, we used F547M/F555W images from seven HST 
pointings that cover the inner ~6 kpc of M83. Our cluster sample was selected from the 
catalogue presented by Silva-Villa et al. (2014), and consists of relatively isolated, well- 
resolved, massive (>10^ Mq) clusters. We use GALFIT to fit the two-dimensional light 
profiles of our clusters assuming an EFF power-law profile. The main results from this study 
are summarized as follows: 

1. The distribution of effective radii is well-represented by a lognormal distribution with 
a median at ^2.5 pc and a range of ~0.3 to 10 pc (Figure 2). Our sample is likely 
incomplete in both very compact and very extended clusters, but the ~3 pc peak is 
present in our data even when including more extended objects. We suggest this is 
strong evidence for a characteristic size of YMGs. Our effective radius distribution is 
very similar to that found by many studies in the literature of both YMGs and GGs 
in the Milky Way and other nearby galaxies. 



2. We find that the mean effective radii, calculated in bins of 0.3 dex in age, tend to 
increase with cluster age (Figure 3). We interpret this as evidence of expansion over 
the hrst few hundred Myrs of evolution. This trend has a power-law slope of 0.26 ± 
0.07. Other studies of young massive clusters in nearby galaxies typically hnd similar 
shallow relationships between age and size, notably in M83 (Bastian et ah 2012), or 
no obvious trend at all. 

3. The effective radii of clusters in our sample do not show a strong correlation with dis¬ 
tance from the center of M83. A weak trend is visible in Figure 4. When we separate 
our clusters into three age bins to remove the age-size correlation in Figure 5, we hnd 
little to no correlation between mean effective radius and galactocentric distance, ex¬ 
cept perhaps for the youngest clusters. Several other studies of young massive clusters 
have found similar results. In studies of GCs, however, some systems display a strong 
correlation between size and distance, such as the Milky Way, and others do not. 

4. We hnd a slight increase in ehective radius with cluster mass, with a power-law slope 
of 0.3 ±0.1 (Figure 6). If we separate the clusters into bins of cluster age, we hnd 
the mass-radius relation becomes stronger with age (Figure 7). Many studies in the 
literature have found little to no correlation between mass and radius for both YMCs 
and GCs. 

5. The distribution of core radii is fairly lognormal in shape with a median at ^1.3 pc and 
a range of ^0.1 to 4 pc (Figure 8). This distribution is less biased against extended 
clusters because of the relaxed rj limit, and it is still strongly peaked. In addition, 
this distribution is somewhat similar to those measured for clusters in nearby galaxies, 
except for perhaps the Magellanic Clouds. We also hnd an overall trend of increasing 
core radius with increasing age (Figure 9), which is very similar to the ehective radius 
trend with age. This trend has a power-law slope of 0.24 ± 0.05. 

6 . We measure the power-law slope of the EFF light prohle, rj, for each cluster. The 
distribution of rj values peaks in the lowest bin and rapidly falls oh to higher values 
(Figure 10). The median rj is 1.5, and the majority of the clusters are best described 
by ?7 < 3.0. Other studies hnd a similar range of rj values as we do here. In an attempt 
to hnd change in rj with other cluster parameters, we divide the rj distribution into bins 
of cluster age, galactocentric distance, and mass (Figure 11). While statistical tests 
suggest signihcant diherences exist between some of the bins in galactocentric distance 
and mass, the diherences are not obvious. Evolution of rj should be further tested with 
cluster population studies in more galaxies. 

7. The expansion of our clusters with age is most likely due to mass loss by stellar evolution 
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or GMC interactions. In addition, the tidal field of M83 appears to play a minor role 
in setting the size of young clusters in our sample. If our clusters were strongly limited 
by tides, we would expect to see a contraction of the effective radius with age instead 
of an expansion. We would also expect to see a stronger mass-radius relation, of order 
Teff oc (Equation 4). We do find a slope of ~0.5 in the 200-300 Myr bin in 

Figure 7, which may be evidence of clusters becoming tidally-limited as they age. This 
is supported by a slight increase of the average rhm/''"j with cluster age. However, in 
conjunction with the lack of a strong relationship between radius and galactocentric 
distance in Figure 5, the evidence suggests that the structure of younger clusters with 
masses >10"^ Mq does not appear to be strongly affected by galactic tides. 
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